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It was shown recently [D.A. Lidar et al., Phys. Rev. Lett. 81, 2594 (1998)] that within the 
framework of the semigroup Markovian master equation, decoherence-free (DF) subspaces exist 
which are stable to first order in time to a perturbation. Here this result is extended to the non- 
Markovian regime and generalized. In particular, it is shown that within both the semigroup and 
the non-Markovian operator sum representation, DF subspaces are stable to all orders in time to a 
symmetry-breaking perturbation. DF subspaces are thus ideal for quantum memory applications. 
For quantum computation, however, the stability result does not extend beyond the first order. 
Thus, to perform robust quantum computation in DF subspaces, they must be supplemented with 



/^ . quantum error correcting codes. 

W^ ! I. INTRODUCTION 

p^ . The power promised by quantum computers H has initiated an intense scrutiny of the physical viability of these 

^ ■ computers [0-0 . The central obstacle in the experimental realization of such computers has proven to be maintaining 
T— I , the quantum coherence of states which form the cornerstone of the speedup promised by quantum computers. The 
^T ' main cause of this degradation of the quantum coherence is the coupling of the computer to the environment, and 
^— ^ , the subsequent decoherence induced by this coupling. To overcome this difficulty, Quantum Error Correcting Codes 
Si^ (QECC) inspired by classical coding theory have been developed ^. These codes are "active" in the sense that 

^^ , decoherence is fought by continuous application of error correction procedures to quantum bits (qubits), which are 
^v • encoded over the Hilbert space of several physical qubits. Another approach has emerged more recently, in which 
■^ I the structure of the physical decoherence process is used to protect the precious quantum coherence [p|-pd| . These 
"passive" error prevention codes rely on symmetries of the decoherence process to encode qubits into states which 
reside in the decoherence-free (DF) subspaces of multiple physical qubit systems. The conditions under which such 
DF subspaces can exist have been established in both the Lindblad (Markovian) formulation [pUlO|] and for the non- 
j^ ' Markovian case H]. 

^ ' The non-Markovian formulation of the reduced dynamics of a subsystem is used extensively in the quantum com- 

O^ putation literature, where it is known as the "operator sum representation" (OSR) |l2J. While it is exact , it is not 
►^ always clear how to separate system from bath in this approach, since for any finite-dimensional bath the dynamics 



are reversible (albeit potentially with a very large Poincare recurrence time) |13 as long as no measurements are 
made. In the Lindblad formulation, on the other hand, the dynamics being described are a priori that of the system 
alone, subject to preservation of complete positivity of the system density matrix. The dynamics are irreversible, but 
the price paid is that the Markovian approximation must be invoked p4| . The first aim of this paper is to clarify the 
relation between these two formulations. We do this with a derivation of the Lindblad semigroup master equation 
(SME) from the OSR in a way which highlights the differences and similarities between the two. Derivations of master 
equations are numerous, starting with the work of Zwanzig and others in the late 1950s pq [: our current derivation is 
novel in that it starts from a fixed-basis form of the OSR which is constructed to formally resemble the SME, yet is still 
exact. This fixed-basis representation allows one to clearly identify the manner in which the temporal coarse-graining 
assumption is invoked in the process of making the transition from non-Markovian to Markovian dynamics, as well 
as the consequences of this. The fixed-basis OSR equation leads to several important results concerning the decay of 
quantum coherence. First, we show that any finite total Hamiltonian will have a zero first order decoherence rate in 
the non-Markovian case, but that this feature can be destroyed by the coarse-grained time averaging made upon going 
to the Markovian limit. Second, a non-vanishing first order decay rate within the non-Markovian OSR formulation 
necessarily implies a singularity of the total Hamiltonian. We show explicitly how this accounts for the behavior seen 
in the prototypical example of phase damping, which is commonly used in the quantum computation literature. 

The second, and main, focus of this paper is the use of the fixed-basis OSR representation to provide a general 
stability analysis of DF subspaces. In Rcf. [|l^ it was shown that, within the SME, DF subspaces are stable to 
symmetry-breaking perturbations, to first order in time. This leads to a lowering of the threshold for fault-tolerant 
quantum computation ||l^,0 in DF subspaces. This stability result was obtained for the "memory-fidelity", i.e.. 



for a quantum computer not subject to external "programming pulses". Since DF subspaces also appear in the 
non-Markovian setting QJl^l j it is natural to enquire then whether the robustness of DF subspaces with respect to 
perturbations extends also to the non-Markovian situation. We make a stability analysis within the OSR here to 
address this, and find that the robustness of the DF subspaces is upheld in the more general situation. Indeed, we 
show that in fact, within both the SME and the OSR, DF subspaces are stable to a symmetry-breaking perturbation 
to all orders of time. 

The generality of this stability result for memory fidelity docs not extend however to the "dynamical fidelity", 
which measures the preservation of quantum coherence in the case of a quantum computer that is subject to external 
programming pulses. We examine this dynamical fidelity here and find that in this case, a DF subspace is stable 
only to first order in time. Thus while a DF subspace can drastically extend the decoherence time for quantum 
memory, operations performed on the DF subspace must be performed rapidly (in comparison to the perturbing 
error rate) in order achieve a similar extension of the decoherence time for implementing quantum computation. 
Barring methods that rely on symmetrizing operations on a time-scale faster than the primary decoherence rate [nSl , 
DF subspaces were, however, never meant to be a complete solution to the problem of decoherence on a quantum 
computer. Their usefulness lies in the elimination of the primary source of decoherence and the subsequent lengthening 
of the decoherence time to one determined solely by the decoherence due to perturbing errors. DF subspaces should 
be supplemented by a QECC in order to achieve a decoherence free quantum computer. This is possible using a 
concatenation of DF and QECC codes, as was shown in Ref. [|l|. The instability of DF subspaces while the system 
is evolving thus sets a lower bound on how rapidly operations on the DF subspace must be performed in order realize 
the robustness of DF subspaces. 



The structure of the paper is as follows. In Sec. [I A a brief review of the OSR formalism is presented, followed by 



the derivation of the fixed basis form for the OSR equation in Sec. II B, This equation is used to derive the SME in 



Sec. [I C. In Sec. Ill we turn to the main subject of the paper, namely short-time expansions of the fidelity. Addressing 
first the fidelity over the entire system Hilbert space, we derive the first order decoherence rate within the OSR and 
show that it vanishes generally in the OSR for non-Markovian dynamics, provided that the total Hamiltonian is of 
a non-singular form. This condition precludes the situation of a system coupled to an infinite number of degrees of 
freedom. We then show that this first order rate may become finite as a result of the coarse-grained time averaging 
performed on going to the Markovian limit. The generality of the non-Markovian result appears initially surprising, 
since there exist elementary examples of non- vanishing first order decoherence rates in the non-Markovian situation. 



In Sec. IIIB, we show with the example of the well-known case of decoherence due to phase damping, how this 
reflects an underlying singularity in the total Hamiltonian. Sec. |^ then deals with the special issue of stability for 
the DF subspaces. After a brief summary of the conditions for DF subspaces in the two approaches, we then show 
that DF subspaces have enhanced stability over the general system Hilbert space, namely that they are stable to all 
orders of symmetry-breaking perturbations, both within the non-Markovian OSR and the Markovian (SME) limit. In 
Sec. M, we then address the "dynamical fidelity" of DF subspaces under external fields, corresponding to a "quantum 
computer program." We show that this is stable to a lesser extent, possessing a vanishing first order decoherence rate 
in both the SME and the OSR, but having non-zero terms of higher order in time. We conclude with a summary and 
discussion of the implications for quantum computation in Section ^. 

II. THE OPERATOR SUM REPRESENTATION 

A. Brief Review 

The dynamics of a quantum system S coupled to a bath B, which together form a closed system, evolves unitarily 
under the combined system-bath Hamiltonian Hsb = H (g) Is -|- I5 ® Hb + H/. Here H, H^ and H/ are, respectively, 
the system, bath and interaction Hamiltonians, and I is the identity operator. Assuming that S and B are initially 
decoupled, so that the total initial density matrix is a tensor product of the system and bath density matrices {p and 
Pb respectively), the system dynamics are described by the reduced density matrix: 

p(0)^^p(t) = TrB[U(/7§5ps)Ut]. (2.1) 

Here Tr^ is the partial trace over the bath and U = cxp(— -IHsst). By using a spectral decomposition for the 
bath, Pb — X^iy '^l^)(^li ^^^ introducing a basis {|n-)}^=i for the iV-dimensional system Hilbert space 7i, this can be 
rewritten in the OSR as ||l|,0: 

K 

p(i) = ^A,(t)p(0)At(t), (2.2) 

i=0 



where the Kraus operators {Ai} have matrix elements given by: 

[^^]mn(t) = ^ {m^V {t)\v)\n) ; i = {^i,v). (2.3) 

Here |i^, /i) are basis elements of the bath Hilbert space, and K ~ Ng, where Nb is the nmnber of bath degrees of 
freedom. Also, by unitarity of U, one derives the normalization condition 

K 



^A|A, = I, (2.4) 

i=0 

which guarantees preservation of the trace of p: 

Tvipit)] - Tr[^ A, p(0) A|] = Tr[p(0) ^ A|A,] = Tr[p(0)]. (2.5) 

i i 

The Kraus operators belong to the Hilbert-Schmidt space A{'H) (itself a Hilbert space) of bounded operators acting 
on the system Hilbert space, and are represented by iV x A^ matrices, just like p. A{'H) is endowed with the scalar 
product 

(A„Aj)=Tr[A,A]]. (2.6) 



B. Fixed-Basis Form of the Operator Sum Representation 



While the OSR evolution equation, Eq. (2.2), is perfectly general, it presents difficulties when trying to separate 
out the unitary evolution of the system from the possibly non-unitary decoherence which occurs from the coupling of 
the system to the bath. The reason is that in general, each Kraus operator will contain a contribution from both the 
unitary and the non- unitary components of the evolution. When one makes the assumption of Markovian dynamics, 
however, we shall see that the semigroup master equation (SME) does separate the evolution of the system into 
unitary and non-unitary parts pO| . This motivates us to manipulate the OSR into a form similar to the SME, but 
without making any Markovian assumption. 

It is convenient for this purpose to introduce a fixed operator basis for A{Ti) |gl[]. Let {Kq,}^^q, with Kq = I, be 
such a basis, so that the expansion of the Kraus operators is given by: 

M 

A,(i) = ^&,o(i)Ko. (2.7) 

For example, {Kq}*^ could be the generators of the Lie algebra su{N) [M = N"^ -- 1 ), or some sub-algebra thereof 
(with M < N^ — 1) |1^ . Under this expansion, the OSR evolution equation, Eq. ( ^.2|) , becomes 

M 
P{t) = Y. Xa/3(0Kap(0)K]3, (2.8) 



where Xaf}{t) is the hermitian matrix 



K 



xcp{t)^Y.^^a{mp{t). (2.9) 

4=0 

Likewise the normalization condition, Eq. ( p.4[ ), is given by. 



M 
Y, xMmlK^ = I. (2.10) 

a, 13=0 



Next we separate out the action of the identity on both Eq. (UM and Eq. (2.10) yielding: 



pit) = XooP(O) 



M 

E 



M 



[Xao(i)K„p(0) + Xoc(i)p(0)Kl]+ ^ XcMOK.pWK^ 



a,/3=l 



(2.11) 



M 



M 



Xool + Y. [xoa(t)Kl + Xao(i)Ko] + Y. XcpimlK^ = I. 



(2.12) 



a,/3=l 



Multiplying Eq. ( 2.12| ) by \p{^) separately from both t he left and the right, adding the resulting equations, and 
substituting the resulting expression for Xoop(O) into Eq. ( 2.11 ), we find: 



M 

p(t)-p(0)--J[S(t),p(0)] + - Y Xa/3(0([Ka,p(0)K^] + [K„p(0),K^]), 



(2.13) 



a,/3=l 



where S(i) is the herniitian operator defined by 



ih 



M 



m^iyYh-o{t)K^-Xoc.{t)Ki]. 



(2.14) 



Q = l 



Eq. ( ^.13| ) is the desired result: it represents a fixed-basis OSR evolution equation. This generally resembles the SME 
in form, but this resemblance should be considered with caution: it can be shown by explicit resummation that for 
purely unitary evolution, p{t) = U(t)p(0)Ut(i), the terms in the fixed-basis OSR are [S(i),p(0)] = [sin(S^) , p(0)] 
and 5Ea/3(---) = U(t)p(0)Ut(i) + (i/n) [sin (Si) ,p(0)] - p(0). Hence the first term alone [S(i),p(0)] does not 
necessarily account for the entire unitary dynamics, as one might naively be led to suspect. 



Comparison of the Fixed-Basis Operator Sum Representation Equation with the Semigroup Master 

Equation 



We recall that in the semigroup approach, under the assumptions of (i) Markovian dynamics, (ii) complete positivity, 
and (iii) initial decoupling between the system and the bath, the system evolves according to the SME ||l4| , |2(]| , P4 : 



dpit) 



L[pit)] 



dt "^^ " h 

M 



[H,p(i)]+Lz,[p(i)] 



^dW)] = \Y. «o/3([Fa,p(i)F^] -f [F„p(i),FJ3]), 



(2.15) 
(2.16) 



a, 13=1 



where 0^/3 is a constant hermitian matrix. This equation bears a clear resemblance to Eq. ( 2.13 ). In fact, taking the 
derivative of Eq. (2.13), we find 



dp[t) 
dt 



-[S(t),p(0)] 



M 

E 

Q,/3=l 



Xaf3 



(t)([K„,p(0)Kt] + [K„p(0),Kl, 



(2.17) 



Noticing the subtle and important differences between the SME Eq. (2.1i) and this OSR evolution equation (2.17) 
allows us to understand the exact manner in which the semigroup evolution can arises from the OSR evolution under 
the above mentioned conditions. The most important difference between these two equations is the fact that the SME 
provides a prescription for determining p{t) at all times t, given p{t') at any other time t' > 0, whereas Eq. ( 2.17 ) 
determines p{t) in terms of p(0), i.e., at the special time t — where the system and the bath are in a product state. 

We now show that explicit use of a coarse-graining over time, together with the above-mentioned assumptions, 
leads one naturally from the OSR evolution equation, Eq. (2.17) to the SME. We note that it is of course possible to 
derive the SME with other methods, (such as adding an infinite bath |g3|,|^) given the appropriate assumptions. Our 
goal here however is not so much to rederive the SME, as to specifically establish a route from the non-Markovian 
OSR to the Markovian SME. 

At this point it is useful to introduce a time-scale r for the bath "memory" (whose definition will be made more 
precise below) and to coarse-grain the evolution of the system in terms of this time scale: 



Pj = pO't); x»i3:j = XQ/aO'r); j e N. 



Further, rewriting the OSR Eq. ( |2.13| ) as p{t) = A(i)p(0) and defining L(t) through A{t) = Texp /^ L(s)ds 
that 



dt 



^L{tMt)]. 



(2.18) 
we find 

(2.19) 



rO + l)Tf 



Defining Lj = J. L{s)ds , with rn ~ t, we have 



L(s)(is 



n-l 

i=o 



(2.20) 



Next we will make the assumption that on the time scale of the bath r, the evolution generators L{t) commute in the 
"average" sense that [Lj,Lfe] — 0,Vj, fc. Physically, we imagine this operation as arising from the "resetting" of the 
bath density operator over the time-scale r. Under this assumption, the evolution of the system is Markovian when 
t »t: 



A{t) - [] exp [tL,] 



Further, under the discretization of the evolution, this product form of the evolution implies that 

Pj+i =exp[rLj] [pj]. 
In the limit of r ^ t we expand this exponential, to find that 



(2.21) 



(2.22) 



(2.23) 



This equation is simply a discretization of Eq. ( 2.19 ) under the assumption that r <C 6>, where 9 is the time-scale of 
change for the system density matrix. Notice in particular tha t the RHS of Eq. ( p. 23 ) contains the average value of 
L(t) over the interval. Now, from the OSR evolution equation (2.17), we know the explicit form of L{t) over the first 
interval from to r. Discretizing over this interval we find that 



pi - Po 



M 



,P0 



- J2 (Xa/3>([K„,poKl5] + [K„p(0),Kt 



a, 13=1 



Lo[po], 



where 



T Jo 



X{s)ds. 



(2.24) 



(2.25) 



Thus, in the sense of the coarse graining above we have arrived at an explicit form for Lq. Consider the evol ution 
beyond this first interval. Deriving an explicit form for Li and for higher terms is now impossible because Eq. ( 2.17 ) 
gives the evolution in terms of p{0). However, since we have made the assumption that the bath "resets" over the 
time-scale r, we expect the bath to interact with the system in the same manner over every r-length coarse-grained 
interval. This is equivalent to assuming that Li = LojVi (which of cour se is the most trivial way of satisfying the 
Markovian evolution condition [Li,Lj] = 0,Vi,j). Then, using Eq. ( 2.23 ), one is led to the well known form of the 
semigroup equation of motion: 



dpjt) 
dt 



M 

J[(s),p(i)] + - ^ (xa;3)([K„,p(i)Kt] + [K„p(t),4 



(2.26) 



Q,/3=l 



(under the natural identification of the K's with the F's of the SME). 
We can write this equation of motion in an alternative form which distinguishes between the system and bath 



contributions to the Liouvillian evolution. Because Eq. (2.17) is linear in the Xctpii) matrix, one can calculate Xagi^) 



m. 



for the isolated system and hence define the new terms which come about from the couphng of the system to the 
bath: 



- ,.(0) 



(1)^ 



x.p{t)^x':i,(t)+x'o.'0{t) 



(2.27) 



The terms which corr espon d to the isolated system will therefore produce a normal — (i/?i)[H, p{ty\ Liouville term in 
Eq. ( |2.26D . Thus Eq. (|2.26| ) can be rewritten as 



dp{t) 
dt 



H 



1 *^ 
S(^^) ,p(i)] + 2 E (^S) ([Ka,p(t)4] + [K„p(t),4]) , 



(2.28) 



Q,/3 = l 



which with the identification of {xap) with aap, and Kq, with F^, is equivalent to Eqs.(2.15)-(2.16), except for the 
presence of the second term derived from ( S ) in the Liouvillian. This second term inducing unitary dynamics on 



the system, ( S*^^^ V is referred to as the Lamb shift. It explicitly describes the effect the bath has on the unitary part 

of th e sys tem dynamics and "renormalizes" the system Hamiltonian. It is often implicitly assumed to be present in 
Eq. (|T5l) H. 

In summary, we have shown in this Section how coarse-graining the evolution over the bath time-scale r allows 
one to understand the connection between the OSR and the semigroup evolution. Specifically, we have made the 
assumptions that (i) the time-scale for the evolution of the system density matrix is much larger than the time-scale 
for the resetting of the bath (r 3> 0), (ii) the evolution of the system should be Markovian ([Li,Lj] = 0,\/i,j), and (iii) 
the bath resets to the same state so that the system evolution is the same over every coarse graining (L^ = Lo,Vi). 
This last assumption can be relaxed, and replaced by a n en semble average taken over the different states to which 
the bath resets, i.e. L^ — ((Lj))- The importance of Eq. ( ^.17 ) lies in the fact that it allows one to pinpoint the exact 
point at which the assumption of Markovian dynamics are made and further, due to the general likeness of its form 
to the SME, provides an easily translatable connection when going from the non- Markovian OSR to the Markovian 
SME. Notice also that the assumption of Markovian dynamics introduces an arrow of time in the evolution of the 
system through the ordering of the environmental states: the system evolves through time in the direction of each 
successive resetting of the bath. 



III. SHORT-TIME EXPANSIONS OF THE MEMORY FIDELITY IN THE OSR 



The mixed-state m,em,ory fidelity [G6J 



F^{t)^Tr[p{0)p{t)] 



(3.1) 



is a good measure of the degree to which a system serves as a perfect quantum memory. Fj^ (t) is the mixed state 
analog of the survival probability for a pure state wavefunction. When the initial preparation is pure, a perfect, 
noiseless quantum memory will have -F'm(i) = 1, but in the noisy case < i^m(^) < 1- If one starts out in a mixed 
state then F,„(0) < 1, and it is usually necessary to resort to some kind of purification |2j. We will consider here 
only short-time expansions of the fidelity, since it is known that using QECC it is possible to restore the coherence of 
a quantum system as long as corrections are applied sufficiently frequently [|l6[ . Thus, we perform a power expansion 
of the fidelity in time p8| 



^-W = E ;!i 



(3.2) 



where the decoherence rates are defined as 

-^{Tr[p(0)p(")(0)]} 

and p*^") denotes the n^^ time derivative of the density matrix. 



l/n 



(3.3) 



A. First Order Decoherence Rate in the OSR 



Throughout the Hterature on decoherence there abound many examples of non-zero first order decoherence rates 
(e.g. ||9|, |29| , p0| ) . Specific attention has been given to maximizing this time scale in order to maintain long-lived coherent 
states. We therefore pose the question here, how do the first order decoherence rates for non-Markovian evolution 
behave within the OSR? The first order decoherence rate is given by 



f 



= Tr[p(0)p(0)]. 



We note that by substituting in the reduced density matrix, Eq. (2.f ), and evaluating the derivative at t = 
immediately led to the vanishing of the first order decoherence rate from the cyclic property of the trace: 



(3.4) 



0, we are 



f 



= Tys[p{0)Ti-bHB.sbp{0) + ip{0)B.sB]] = 0. 



(3.5) 



Therefore for the general non-Markovian dynamics, first order decoherence rates are rigorously zero, provided that 
tlsB is finite. What is not obvious from this simple manipulation is how a coarse-graining procedure can lead to the 
commonly encounte red n on- vanishing first order decoherence rates. This therefore suggests employing the "pre-coarse- 
grained" OSR, Eq. ( 2.f7 ), fo r th e derivative, and then carrying out the specific coarse-graining procedure outlined in 
the previous sectio n on Eq. (3^), in order to understand how non- vanishing first order decoherence rates can arise. 
Using Eq. (2.17) the first order decoherence rate becomes 



1 
— =Tr 

n 



/ . M 

p(0) -J[S(0),p(0)] + - Y. Xa/?(0)([K„,p(0)4] + [K„p(0),4] 



a, 13=1 



(3.6) 



Using the decomposition of the Kraus operators, Eq. (2.2), and knowing that U(0) = I5 (g) Ig, we find that Ai(0) = 
V^Is<5i,(i/.iy)- Thus, since the Kq,'s form a linearly independent basis, it follows, using Eq. (p/7|), that the expansion 
coefficients must be 



biaiO) = SaoVvS, 



ijlu^iy) • 



By direct evaluation, 



Xa/3(0) = X! ^ (^aob*„^^),piO) + '5/30&(i.,^),a(0) 



(3.7) 



(3.8) 



which implies the vanishing [as long as 6(,/^iy),Q(0) remains finite] in Eq. (3^) of every term except Tr[p(0)[S(0),p(0)]]. 
However, this in turn vanishes by cyclic permutation of the trace. Thus we see that within the OSR, the first order 
decoherence rate is always zero when the &(i/,i.),a(0) remain finite. To determine the significance of this restriction we 
choose as a basis for the Kraus operators a set of M Kq, 's which form a Lie algebra and hence have a suitably defined 
inner product [Eq. (p!q)]: 



Tr 



KqK^ 



_ / (5„^ for a>l,(3>l 
NSp^Sao otherwise. 



We then find using the definition of the Kraus operators, Eq. (2.3 ), 

\..)At) = Tr [^iMv\\J{t)\v)] . 
Differentiating this and recalling that U(t) — exp(— iHssi/?!), we find that 



(3.9) 



(3.10) 



&(..),a(0) = -^Tr [K.i^{v\iisBmy) 



(3.11) 



Thus in order for bi^^\,^{Q) to remain finite, (z^|H5B(0)|i^) must be finite. Hence our conclusion that the first order 
decoherence rate vanishes in the OSR is valid for any finite total Hamiltonian (by which we mean its matrix elements), 
and conversely, any finite total Hamiltonian will have zero first order decoherence rate pll . 



Examination of our derivation of the SME, Eqs. ( 2.26| ) and ( ^.281 ), now shows how non-zero first order decoherence 
rates can arise when the evolution is considered to be Markovian. In the derivation of the semigroup equation in 
the Markovian hmit we made the as sump tion that the matrices Xapit) can be identified with the constant mat rices 
ttap of the semigroup equation, Eq. ( 2.16 ). However, when this i s do ne, the matrix elements Xap{0) in Eq. (B.C) are 
replaced by their time-averaged values, for which the relation (3.8) no longer applies. Hence, in general, the first 
order decoherence rates are necessarily not zero when the Markovian coarse-graining is applied. For a finite total 
Hamiltonian Hsb, non-zero first order rates are therefore seen to be an artifact of the Markovian assumption, and 
their appearance emphasizes the delicate nature of the transition to the Markovian regime. 



B. Example: Phase Damping 

The restriction to a finite total Hamiltonian above may at first sight seem obvious. However, consider, for instance, 
the often quoted example of phase damping of a qubit. In this case, it would appear that there is a finite first order 
decoherence rate. Yet, it is often presented within the OSR ||3^,y , which, as we have just shown above, would predict 
zero first order decoherence rates for any non-singular Hamiltonian. In this example, the Kraus operators are given 

by[il 



An = 



1 




-At 



Ai = 









Vi~^ 



-2Xt 



(3.12) 



and a simple calculation using these operators yields a minimum first order decoherence rate of 1/ti = —A/2. How 
can this be? 

To resolve this dichotomy, we consider how the above phase damping Kraus operators are generated from the 
unitary dynamics of a qubit system S and a qubit bath B. The evolution operator 



V{t) 






V -VT 







-2Xt 







VT 



-2Xt 



\ 



-Xt 



J 



iO) 

TO) 
Ti) 



(3.13) 



[where the first qubit represents the bath (| t), 1 1)) and the second represents the system (|0), |1)) as denoted in the 
columns above] with the bath initially in the state 1 1), immediately gives the Kraus operators of Eq. ( 3.12| ). Now, it 



is easy to calculate the Hamiltonian which generates this evolution, [using llsBit) — ih J^ Ut(t)] 



Us Bit) 





~g{t) 



g{t) 



(3.14) 



where 



g{t) = ih 



76 



-7* 



VT 



-27* 



(3.15) 



However, we see that as i — > 0, g(t) — > cxd. Thus, in this simple example, we find that at t = 0, the Hamiltonian 
becomes singular. This illustrates our claim that first order decoherence rates in the OSR are the result of an infinite 
Hamiltonian, and do not contradict the general OSR result of zero rates for finite Hamiltonians. 

The diverging Hamiltonian in this example is in fact equivalent to non-closedness of the system S + B. It is 
well known that phase damping can be generated by a model of random phase- kicks [ p4[ , which implies an external 
random force, i.e., that the system S -\- B is in fact not closed. Since this is in contradiction to our initial assumptions 



(Sec. [I A), it should not come as a surprise that a non- vanishing first order decoherence rate is found in this situation. 
A similar divergence will result of course from a bath with an infinite number of degrees of freedom [g3|,|^. The OSR 
phase damping example thus can still be used (as is commonly done in the analysis of quantum error correction) 
under the caveat that one cannot claim that it arises from a finite closed system (S+B). 



IV. EFFECT OF SYMMETRY BREAKING PERTURBATIONS ON MEMORY FIDELITY OF 

DECOHERENCE-FREE SUBSPACES 

Our discussion in the previous sections was completely general, dealing with the decoherence of the entire system 
Hilbert space. We now restrict our attention to the behavior of the fidelity in DF subspaces. We first briefly 
summarize the basic theory of DF subspaces and then generalize the first order stability results obtained within the 
SME in Ref. |§. 

A. Theory of Decoherence-Free Subspaces: Markovian vs non-Markovian Approach 

Recently, conditions for the existence of decoherence free subspaces within the framework of the Markovian SME 
approach P,pi| and in a non-Markovian [R setting were derived. We first clarify here the connection between the 
SME and the non-Markovian results. 

In the SME approach it was shown that a necessary and sufhcient condition for decoherence free dynamics 
(L/)[p(i)] = 0) in a subspace 7i = Span[{|i)jJ'-^}] is that all of the basis states |«) satisfy the condition 

Fa|?)=c„|?) Va,? (4.1) 



where the Fq,'s are the error generators in the semigroup Eq. ( 2.16| ). Since the {Fq,} form a Lie algebra £, this 



condition has a simple group-theoretic interpretation, namely, the DF states are th e sin glets of £, i.e., they a re th e 
states that transform according to the 1-dimensional representations of iZ . From Sec. [I C, and in particular Eq. ( ^.26 ) 



it follows that these error generators become identical to the {Kq,} (the basis operators in the fixed expansion of the 
Kraus operators) when the short-time averaging approximations leading from the OSR to the SME are made. 

Within the framework of non-Markovian evolution, it has likewise been shown 0,^ that a necessary and sufficient 
condition for decoherence free dynamics over a similar subspace H = Span[{|i)^J'-^}] is that all of the basis states \i) 
satisfy the condition 

Sa\t)=Ca\t) yaj (4.2) 

where Sq's (system operators) are defined by the interaction Hamiltonian 

H7 = ^S„(^B„. (4.3) 



The Bq are bath operators. The fixed-basis OSR equation (2.13) sheds light on the relationship between these two 
DF conditions. To the extent that the error generators Fq can be derived from expansion of the Kraus operators 
with a subsequent short-time averaging approximation, the DF condition given for the non-Markovian dynamics are 
more general than that given by the semigroup approach. We notice that in the limit of small averaging time r, the 
error generators F^'s will, in fact, correspond directly to the Sq's. This can be seen by expanding the full evolution 
operator U(t) to first order in r: 



'— -'2^ 



U(r)=I--HsBT + 0(T^). (4.4) 

To this order, the Kraus operators Ai{t) [see Eq. ( pT^ )] will only contain terms which correspond to terms that 
appear in the Hamiltonian Hsb- These are H, the system Hamiltonian; I, the identity; and the Sa's from Eq. ( [4.3| ). 
Terms corresponding to H will result in unitary evolution of the system, while I is removed from the set of Kraus 
operators (since these are considered error generators) via our derivation of the SME. To first order in time, therefore 
the only error generators are the Sq's. Thus in the case of small averaging time r we see that the two DF conditions 
are exactly equivalent. 

However, it is important to note that the SME approach has other advantages. Thus in many cas es it is either 



impractical or undesirable to derive the {Fq.} from a short-time expansion of the type discussed in Sec. II C . In fact, 
in Lindblad's axiomatic approach uM, the {Fq} are the primary objects and they do not follow from an e xpa nsion 



of a unitary operator. While the {Fq} are often identified heuristically from a factorization such as in Eq. (O) [^ 
in some cases (notably strong coupling) one simply cannot clearly separate system and bath in the form assumed in 
that equation. In this sense, then, the SME provides greater generality than the non-Markovian approach within the 



Hamiltonian representation subject to Eq. (4.3). Motivated by this aspect, we take condition (4.1) to be necessary 
and sufficient for DF subspaces. 



Further, one should note that while the semigroup DF condition, Eq. (4.1), guarantees that the evolution of the 
system will be unitary, the system may still be subject to unitary evolution induced by the bath in the form of the 
Lamb shift. Such bath-induced evolution, although it does not introduce decoherence, is undesirable in the course of a 
quantum computation. To the extent that we desire the DF subspace to serve as the basis for a quantum computer, we 
therefore must impose one of two conditions on the DF subspace. These conditions are: (1) suitable control over the 
system is obtained so that the Lamb shift term can be canceled out, or, (2) the Lamb shift docs not induce dynamics 
on the subspace, i.e., HLamb|«) — hi\l). Under the first condition, the DF subspace is not reduced in dimension. 
However this condition may be physically impossible to realize. The second condition does not make any assumptions 
about the amount of control one has over the system Hamiltonian, but it may have the undesirable effect of causing 
a reduction in the size of the DF subspace. Similar conclusions regarding the effect of the Lamb shift hold for the 
non-Mar kovian case. 

Finally, in both the SME representation and the OSR, the implicit assumption has been made that the system 
Hamiltonian does not induce the evolution of states from within DF subspace to states outside of the DF subspace: 



mi) 



E'^^^i^)- 



(4.5) 



It is a simple matter to satisfy this additional condition in the SME. However, this is not the case in the OSR due to 
the absence of an explicit appearance of H. 

We now proceed to the question of stability of DF subspaces in the non-Mar kovian and Markovian cases. The 
calculations in the next two subsections are rather tedious, and the reader who is not interested in the details may 
wish to skip directly to the results for the memory fidelity, summarized in Table I. 

B. Stability of the Memory Fidelity in the Non-Markovian Case 

Consider the addition to a DF subspace of new perturbing terms in the interaction Hamiltonian: H'^^ = Hgs+eHj. 
In this case, we find that to first order in e , the new full evolution operator is given by 



u'w = E 



{-it/ny 



(H 



SB 



eWr 



n=0 



U(t) 



k—l n^k 



ti^/,'«(H..,H;) 



(4.6) 



where 



n 



(1 



/: 



(2 



(Hss,Hj) = H'j 

(Hss, Hj) — HssHj + HjHgB 

(HsB, Hj) = HggHj + HsbHiH-sb + Hj^sb 

(Hss,Hj) = Hj 



(Hss, Hj) = HssHj + HjHsbUi 



H/ Hsb, 



(4.7) 



etc. Here U(t) is the unperturbed evolution operator. From Eqs. (2^) and ( p.7| ) we thus see that to first order in e 
the operators Kq in which the Kraus operators are expanded, will have new terms due to {/A }, hereby denoted by 
{eGp} ^-^, which are proportional to e. These terms modify the evolution over the DF subspace [Eq. ( ^.13 )] so that 
dp/dtt^ dp'/dt = dp/dt + L'(i)[p(0)], where: 



M P 

L'itMO)] ^ - ^[S'{t),pm + ;j EE^"pWLk.,.g,[p(0)] +x;p(i)L, 



h 



G„,K£ 



[pm 



Here: 



a— 1 p— 1 



ih y—^ 



S'it)^e-J2[Xpoit)Gp~xop{t)Gl] 
p=i 



(4.8) 



(4.9) 
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and: 



Lx,y[p] = [x,py'f] + [xp,y^]. 



(4.10) 



Terms of 0(e), not written out explicitly in Eq. (4 



include LeGp,cGjp(0)], U^a^.KjpiO)], etc. Assuming that e < 
1, we may neglect these terms. Now, for the purposes of argument, we will assume that the system has perfect quantum 
memory over the DF subspace in the absence of the perturbing error generators, i.e., F^a (t) = Tr[p(0)/5(t)] = 1. The 
perturbation, however, decreases the fidelity below this perfect value. The modified memory fidelity can be written 
formally as: 



A;— 1 n— fc 



KT' 



Jk) 



where r„ ' represents the 0(e) contribution to the n order decoherence time r„. It was shown in Ref. [hO| that, 

within the SME, the term which is first order in both the perturbation and time vanishes: 1/t} — 0. This left 
open the possibility of terms of order O(et^) and higher spoiling the fidelity. Here we will generalize this result in the 

non-Mar kovian case [the Markovian case will be dealt with in Sec. ( p^V C| )] and show that in fact l/ri = for all 
n, so that the entire 0(e) contribution vanishes, and only terms of order 0{e^t^) can spoil the memory fidelity. For 

simplicity of notation, since we are considering here only the 1/t„ decoherence rates, we drop the (fc) superscript 
from now on. 

The perturbed decoherence rates are thus given by: 



1 



,(») 



(n-l). 



- = Tr[p(o) {p'{t)r'] = Tr[p(o) {l' mmmio 



(4.11) 



Using Eq. (4.8) and noting that the terms involving S' vanish directly by permutation under the trace, we obtain to 
first order in e: 



1 



M P 



EE^"7'ni)Tr[p(0)LK„,.G,[p(0)]]+x:i"-'Ht)Tr[p(0)L,G„Kj/5(0)]] 



(4.12) 



To evaluate this, we need to know Kqp(O). Now, when we expand the Kraus operators about a fixed basis Kq as in 
Eq. (2.7), this basis will consist of all possible products of the three terms H, Sq, and I [recall the definition of the 
Kraus operators, Eq. ( |2.3| )]. Assuming a perfect quantum memory 

(4.13) 

Ca\i) , for a given product of S^'s and H's we can 



[H,p(0)]=0, 

we can commute H with /5(0) and, using the DF condition Sq|?) 
replace each Sq with its eigenvalue c^. Thus, for example. 



H^SiHp = H^SipH = H^cipH = ciH^p. 
It follows that for a DF subspace in the OSR, the basis operators K^ will satisfy the condition 

K„p(i) = d„H'""p(i), 



(4.14) 



(4.15) 



with rria an integer and da a real number. Using this result, we then cycle the trace in Eq. (4.12), and again using 
[H, p(0)] = 0, we find for the first trace 



Tr[p(0)LK„,eG,[p(0)]] 

= ed„Tr[2p(0)H""/5(0)Gt - /5(0)GtH™°p(0) - p(0)Gtp(0)H'^ 
= 0. 



(4.16) 



The second trace vanishes similarly. Thus we see that all orders of decoherence rates must vanish to first order in 
e, i.e., l/r„ — 0. Examining Eq. (4.8), it is clear that the second order e^ term does not lead to a similar vanishing 
of the traces. Therefore we have proved that: under non-Markovian evolution, DF subspaces are completely stable to 
first order under a symmetry breaking perturbation, where by "completely" we mean explicitly stable to all orders of 
time. 
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C. stability of the Memory Fidelity in the Markovian Case 



The stability of DF subspaces with respect to symmetry breakingperturbations in the non-Markovian case derived 
above is an significant extension of the stabiUty derived in Ref. |1^ for the Markovian SME. However, the resuh 
presented in Ref. |1^ only examined the effect of a symmetry breaking perturbation on the first order decoherence 

rate (1/r} — 0). Here we show that the stronger non-Markovian result derived above (1/tA = Vri > 1) also holds 
in the Markovian SME. 

The effect of perturbing a DF subspace in the SME is the addition of new error generators, hereby denoted by 
{eGp} ^^, to the master equation (which was partially treated in Ref. pC|) as well as a perturbing Lamb shift term 
in the master equation (which was not treated in Ref. [Q). These terms modify the SME [Eqs. ( ^.15 )-( 2.16| )] so that 
9p(t)/c>i = L'[p(t)], with 



L'[p(i)]^L[p(i)]--[6HUb,PW] 



P M 



J E E(3P"L^G,,Fjp(t)] + 5;aLF„,eG,[p(i)]) + 0(e2 



(4.17) 



P— 1 Q— 1 



where L[p(t)] is the u nperturbed SME term given by Eq. ( ^.15|) , (i^\^^^\^ is the perturbing Lamb shift, and Lx,y[/9(i)] 
is given by Eq. (4.10). The perturbed decoherence rates are given by 



Tr[p(0)p(")(0)]^Tr[p(0){(L'r[p(0)]}] 



(4.18) 



where (L')"[/5(0)] = L'[L'[...L'[/5(0)]]], n times. To evaluate this expression, recall (i) the DF condition L[/5(0)] = and 
(ii) that we are working only to first order in e. Now, for simplicity, consider first (L')^[p(0)], and denote the second 
and third terms on the RHS of Eq. ( |4.17l ) by Ai[p(i)] and k2\p[t)\. By the DF condition, L[L[p(0)]] =Ai[L[/5(0)]] = 
A2[L[p(0)]] = 0. Also, A4A-,[/5(0)]] is of 0(e ). Only the two terms with L acting on Ai[p(0)] and A2[p(0)] do not vanish. 
This reasoning generalizes easily for n > 2, so we find that to first order in e. 



Tr 



p(0)L 



ra-l 



P M 



eHLmb,P(0)] + i^f](gp„L,G,,F„[/5(0)] +5;aLF„,eG,[p(0)]) 



p^l a—1 



(4.19) 



Now, the DF condition, Eq. (HjJ), implies that the DF error generato rs co mmute with the DF density matrix: 
[Fc,p(0)] — 0. We also again assume perfect quantum memory [Eq. ( 4.13| )]. Thus, for an arbitrary operator 
A G^(7i), we find that we can commute the initial density matrix through the operator L : 



so that: 



L[p(0)A] = - J[H,p(0)A] + IJ2 «c./3(2F„p(0)AFl3 - F;jF„p(0)A-p(0)AF]3F„) 

M 

= -^p(0)[H,A] + - J2 aa/3/5(0)(2F„AF|j-F|5F„A-AF)3F„) 



a,/3=l 



= p(0)L[A], 



l("-i)[/5(0)A] = p(0)l("^i)[A]; l("-i)[A/5(0)] = l("-i)[A]p(0). 

Examining the contribution to the decoherence rates from the Lamb shift term, we thus find that 



1 

'n ' Lamb 



:Tr[p(0)[L("-i)[eHLa.„b],p(0)]]=0 



(4.20) 



(4.21) 



which vanishes by cycling under the trace. Next, the contribution to the decoherence rates due to the symmetry 
breaking perturbing error generators in Eq. ( 1.19|) are given by 
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P M 



(1) ; 2 

Tn / SBP ^ p=l Q=l 



Expanding the first of these terms, 

P A/ 



^ E E Tr {p(0)L('-i)[g,„L,G,.Fj/5(0)] + g;„LF„..G,[/5(0)]]} 



2 E E 5p"Tr [p(0)l("-i) [2Gpp(0)Fl - FlGpp(O) - p(0)Ft.G, 

P A/ 

I E E 5poTr [p(0) (l("-i) [2GpFl]p(0) - l("-i) [FlGp]p(0) - p(0)l("-i) [F^G^] 



= 1 Q = l 

P A/ 



^ E E 5paTr [p(0)l("-i) [[G„ Fl]]p(0) 

p^l a— 1 
P M 

^E E 5p«Tr [l("-i)[p(0)[Gp, Ft,]p(0)] 

p—1 a— 1 



(4.22) 



(4.23) 
(4.24) 



Using F]jP(0) = c*p(0), we see that p(0)[Gp,Fjj]p(0) = p(0)(GpC* - c*Gp)p(0) and thus this term vanish. Similar 
reasoning implies the vanishing of the second term in Eq. (4.22). Thus we have proven that 1/t„ = 0: under 



Markovian evolution, DF suhspaces are completely stable to first order under a symmetry breaking perturbation |36 



V. THE DYNAMICAL FIDELITY 

The results derived in the previous Section imply that DF subspaces are robust to small perturbations when the 
DF subspace is operating as a quantum memory. In order to address what happens when perturbations are made 
on the system as it evolves according to some desired quantum computation, we have to first define an analog of the 
mixed-state memory fidelity for an evolving system. This is 



F^it) = Tv[pu{t)p{t)], 



where pu{t) is the desired unitary evolution 



Pu{t) = Us(i)p(0)U^(i), with Us(i) = exp 



-Hct 



(5.1) 



(5.2) 



Here H5 is the system Hamiltonian, and may include a "program" Hamiltonian which implements a quantum algo- 
rithm on the system. This dynamical fidelity is a good measure of the difference between the desired evolution of the 
system and the actual, "noisy" evolution. Thus, < Fd{t) < 1, with F^{t) = 1 if and only if the evolution is perfect, 
i.e., p{t) = pu{t). 

The decoherence rates for the dynamical fidelity are defined in the same manner as for the memory fidelity: 



"^m-T.^ (£}"-■ ^ = {M«,«wo)'"'i} 



lA 



(5.3) 



A. Markovian Case 



First we consider the dynamical fidelity within the context of the Markovian limit, using the SME approach. We 
restrict our attention as before to DF subspaces, so the density matrix p satisfies Eq. (2.15) with LD[/5(i)] = 0. We 
then imagine this DF subspace to be perturbed by a symmetry-breaking perturbation: L/j t-^ L'jy, where the perturbed 
density matrix satisfies the following SME: 



dp 
dt 



[H5,/5(t)]+L^[p(t)]. 
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Similarly to Eq. ( p~q ) (see also Ref. ||l^), the new terms in this SME are given by: 

M P 



LD[/5(i)] = E E KpLF„,.G,[/5(i)] + a:pL,G„F„[/5(i)]) + 0{e^). (5.4) 



a— 1 p— 1 

The perturbed first order dynamical fidelity decoherence rate is given by 

1 



= Tr 






(5.5) 



The first of these terms vanishes via the arguments given for the memory fidelity (essentially, since F^p — Cap by 
the DP subspace property). The second term also vanishes, by permutation of the trace after using dpu{0)/dt = 
—ih[Hs,p{0)]. Thus we find that DF subspaces are stable to first order in time also when the system is allowed to 
evolve. 

Further, it is easy to see that the higher order dynamic fidelities now do not vanish. For example, the second order 
dynamic decoherence rate contains terms like [H5, p{0)]L'^[p(t)]t=o, which do not allow the simple permutation of the 
trace. 

B. Non-Markovian Case 

Is th ere an analogous result for the non-Markovian situation? We can address this within the OSR. The arguments 



of Sec. Ill A showed that the first order decoherence rate will always vanish, and so according to the arguments given 



above, this applies also to the dynamical fidelity. Examining the second order decoherence rates, we find: 
^^ ~ ^r^2[H5,[H5,/5(0)]]p(0) + [Hs,p(0)]L'(t)[p(0)]*.o + /5(0)|^[/5(0)]| 



Tr 

T2 



(5.6) 



where L' is now given by Eq. (4.8). As in the semigroup analysis above, the second of these traces does not vanish 
(nor does the first, but it will be cancelled due to a contribution from the third). 

Thus we find that for the dynamical fidelity, the effect of a symmetry breaking perturbation results in a second 
order time instability in the system for both Markovian and non-Markovian situations. It is remarkable that while the 
memory fidelity is completely robust to first order (e) perturbing errors, the dynamical fidelity does not show the same 
robustness, with instability arising at second order in time. In terms of real quantum computation therefore, this 
implies that DF subspaces must be supplemented by QECC in order to be truly useful beyond merely providing high 
quality quantum memory. In particular, the instability in the dynamical fidelity implies that in order to realize the 
robustness of DF subspaces to symmetry breaking perturbations, operations on the DF subspace must be performed 
over a time scale short in comparison with the perturbing error rate. Thus, for the realistic scheme in which DF 
subspaces are supplemented by quantum error correcting codes (by, for example, concatenating the DF subspaces 
within QECC as in Ref. ||ll|), if the operations performed on the DF subspace in order to execute the QECC are 
executed frequently and rapidly, the full scheme can provide a significant improvement over pure DF subspaces. 

VI. SUMMARY AND CONCLUSIONS 

We have shown here how the formally exact operator sum representation (OSR) for the time evolution of the density 
matrix can be cast in a form which bears a significant resemblance to the semigroup master equation (SME), through 
the introduction of a time-independent (fixed) operator basis. By using this fixed-basis OSR equation, we were able 
to easily calculate the fidelity in the OSR, as well as to provide a derivation of the SME which makes explicit the 
role played by the coarse-graining assumption. Somewhat surprisingly, we found an important difference between the 
OSR and the SME, namely, the first order decoherence rate always vanishes in the former, but not always in the 
latter, for a finite Hamiltonian. This effect is readily traced to the coarse-graining time averaging assumption within 
our derivation. This result is significant for both error-correction schemes aimed at improving the fidelity, and for 
commonly used simplified models of decoherence. We illustrated the latter with the well-known case of phase-damping. 

Using the fixed-basis OSR, we have then undertaken a detailed study of short-time expansions of the mixed- 
state fidelity under a variety of conditions. The mixed-state fidelity provides a measure of the extent of quantum 
coherence in the system. We have examined both the usual mixed-state memory fidelity relevant to quantum memory, 
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and a dynamical fidelity which we defined to act as a measure of coherence for a time-evolving system. Our main 
achievement here has been to extend the robustness results of Ref. |1^ regarding decoherence-free (DF) subspaces. 
For the preservation of quantum memory, we showed that in both the OSR and SME approaches, DF subspaces are 
stable to all orders in time to a symmetry-breaking perturbation. The first errors entering the quantum memory 
can therefore only be O(e^), where e measures the strength of the perturbation. This result goes beyond the first 
order stability result of 0{et^) arrived at in Ref. [|lO| , which was restricted to the Markovian case. It shows that DF 
subspaces are indeed ideal for quantum memory in all situations. In making this statement, we note that we have not 
shown here how to perforin input and output to the DF subspaces. Further investigation is needed to address this 
issue. The quantum memory stability results are summarized in Table m For the dynamical fidelity, a weaker result 
is obtained, namely that this has only a vanishing first order decoherence rate under a perturbation. Thus the first 
errors entering the dynamical fidelity can be 0{et^). 

This stability analysis of static and dynamic fidelities is of particular relevance for practical implementations of 
quantum computation. The complete stability to perturbations of static fidelity within DF subspaces is very encour- 
aging for use in quantum memory. Thus passive error correction appears to be sufficient for this. In contrast, the 
weaker first order stability condition derived for the dynamical fidelity within DF subspaces implies that application of 
active error-correction techniques will likely be necessary in order to preserve coherence during computation. Further, 
if the operations performed on the DF subspace in order to execute active error-correction are executed rapidly, then 
the full DF stability to symmetry breaking perturbations can be recovered. Such a scheme of combining active and 
passive error-correction is possible by concatenating codes constructed from DF subspaces with active quantum error 
correction, as demonstrated in Ref. |ll[ . 
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TABLE I. DiflFerent order decoherence rates for the mixed-state memory (r) and dynamical (r) fidelities under the various 
conditions cited in the text. Unless otherwise noted, n > 1. 

SME OSR 

General l/n / l/n = 0=" 

l/r„ / 0, n > 2 l/r„ / 0, n > 2 
DF subspaces l/r„ = 0*^ l/r„ — 

memory fidelity for e-perturbed DF subspaces 1/t„ — 1/t„ — 

dynamical fidelity for e-perturbed DF subspaces l/n — l/n — 

l/f„ / 0, n > 2 l/f„ 7^ 0, n > 2 

'^finite bath and finite total Hamiltonian only, 
both with and without perturbing Lamb shift. 
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